Bose-Einstein Condensates in a Cavity-mediated Triple-well 
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We investigate the energy structures and the dynamics of a Bose-Einstein condensates (BEC) in 
a triple-well potential coupled a high finesse optical cavity within a mean field approach. Due to the 
intrinsic atom-cavity field nonlinearity, several interesting phenomena arise which are the focuses 
of this work. For the energy structure, the bistability appears in the energy levels due to this 
atoms-cavity field nonlinearity, and the same phenomena can be found in the intra-cavity photons 
number. With an increase of the pump-cavity detunings, the higher and lower energy levels show a 
loop structure due to this cavity-mediated effects. In the dynamical process, an extensive numerical 
simulation of localization of the BECs for atoms initially trapped in one-, two-, and three-wells are 
performed for the symmetric and asymmetric cases in detail. It is shown that the the transition 
from oscillation to the localization can be modified by the cavity-mediated potential, which will 
enlarge the regions of oscillation. With the increasing of the atomic interaction, the oscillation is 
blocked and the localization emerges. The condensates atoms can be trapped either in one-, two-, or 
in three wells eventually where they are initially uploaded for certain parameters. In particular, we 
find that the transition from the oscillation to the localization is accompanied with some irregular 
regime where tunneling dynamics is dominated by chaos for this cavity-mediated system. 

PACS numbers: 42.50.Pq, 03.75.Lm, 05.45.-a 



I. INTRODUCTION 



Cavity-mediated Bose-Einstein condensates system 
(BECs) has been widely investigated as a promising plat- 
form to explore the exotic many-body phenomena in re- 
cent years both in theory and experiment due to the con- 
trollability of both the confining potential geometry and 
the interatomic interaction [H, Hj. If the cavity resonance 
and the atomic resonance are in the large-detuning limit, 
the dispersive regime can be realized. At this stage, the 
strong coupling between the atoms and the cavity field 
will induces an additional optical potential for the atoms 
and exerts significant mechanical forces on the motion 
of the atoms; meanwhile, the atoms imprint a position 
dependence phase shift on the cavity field vice versa. As 
the field adjusts accordingly, the field causes a back ac- 
tion on the atoms, resulting in the nonlinear atom-field 
interplay [H-IH . This highly nonlocal nonlinearity is quite 
different from the usual local atom-atom interactions. As 
the cavity-atom interaction induces an optical lattice for 
the atoms, then the cavity-mediated optical lattice model 
is one of the usual platform for discussing the cavity- 
mediated ultracold atoms and BECs ■ Based on this 
nonlinear interactions, the cavity field coupling either to 
the center-of-mass pd] - [l3| or to the spin degrees of free- 
dom [l4[ and to the both cases [HI, [l6| i.e., the non- 
linear coupling among the external and internal states 
arc demonstrated. Recently, rapid proposals have been 
intensely investigated, which gives rise to a number of 
interes ting nonclassical features such as optical bistabil- 
ity [l7l - [22| |. collective dynamics [23| and quantum phase 



transitions [24l.l25|. 

On one hand, a cavity-mediated double-well potential 
is the simplest case of a multiwell potential for studying 
the cavity-mediated BECs system, a plethora of novel 
phenomena, such as strong population transfer, the com- 
plex and diverse stationary behaviors and the pseudo- 
sclf-trapping effect, have been observed and analyzed in 
great detail for the this system because of the unique 
cavity-mediated nonlinearity in the context of cavity 
quantum electrodynamics @, [2(| |27| . The role of dis- 
sipation in the system of a light field coupled with ultra- 
cold atoms trapped in a two-well optical potential in a 
lossy cavity is also investigated [28j . Special asymmetric 
double-well case is studied too, the results show that the 
performed numerical simulations on bistability in optical 
lattice have good agreement with the special asymmetric 
double- well case described by a two- mode model [22j . 

On the other hand, If one wants to make a extension of 
the double-well system, the triple well is the straightfor- 
ward natural one. Compared with the double well sys- 
tem, the difference is essential both in mathematics and 
physics for the triple well system, which has the following 
distinguished characteristics. A triple well potential has a 
much richer structure p9j and can supports the simplest 
model for exhibiting the effect of next-nearest-neighbor 
coupling on quantum tunneling [30j , the quantum tun- 
neling may happen between several wells simultaneously, 
then the tunneling dynamics in the triple well will show 
more interesting behavior [3l|. The study of the triple 
well system will also provide a bridge between the simple 
double well and the multiwell systems, as a protype of op- 
tical lattices, helping us provide a bottom-up understand- 
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ing of mechanisms operating in the infinite optical lattice 
[32| . Moreover, similarly to the analogue of the double 
well to the Josephson junction, the triple well is the min- 
imum system that can model the source-gate-drain junc- 
tion, and draws much attention from the perspective of 
atomtronics. Considering these reasons, recently, then 
there is a growing substantial attention in the dynam- 
ics of ultracold atomic clouds in triple-well potentials, A 
variety of proposals to achieve controllable atomic dy- 
namics based on the triple well have been presented ex- 
tensively by employing the standard models of the Gross- 
Pitaevskii equation and the Bose-Hubbard Hamiltonian 
mainly. It includes the eigenstates and the tunneling dy- 
namics j3^ - l35| . the transistorlike behavior (36[, the adi- 
abatic transport of a hole [H?} , and the mesoscopic quan- 
tum superpositions [30j | . However, unlike the optical lat- 
tice and double well cases, no work has been studied in 
the triple well coupled to an optical cavity so far and a 
fundamental question of whether or not loop structure 
and chaotic exist is still open. Then theoretical predic- 
tions on the dynamics of cavity-mediated BECs system 
in a triple well are highly desirable. On the contrary, as 
we will show below, the cavity-mediated BECs system 
does introduce qualitatively novel physics in triple well 
systems. We will discuss how the nonlinear interaction 
leads to s loop structure in the energy band. In addition, 
we also show how this nonlocal nonlinearity leads to a 
transition quantum dynamics characterized by striking 
new phenomena, including the chaotic quantum dynam- 
ics. 

This paper is organized as follows. In Sec. II the 
mean-field approximation is employed to map the quan- 
tum model into a classical Hamiltonian describing the 
atoms and the corresponding dynamical equations are 
derived. In Sec. Ill we study the optical bistability in 
the atom-cavity system and achieve how the intra-cavity- 
photons and energy band structure depends on the cav- 
ity detunings. In Sec. IV we investigate the transition 
to self-trapping of BECs in one well and in two wells, 
respectively, showing some interesting regions that the 
chaotic motion and steady behavior emerge alternately. 
The topics are extended to a tilted triple-well system in 
Sec. V. Finally, the results are briefly summarized in 
Section VI. 



II. THE BASIC MODEL 



Consider a BECs system confined to in a triple well 
coupled to an external single mode of high-finesse opti- 
cal cavity (schematically sketched as in Fig.(l)). This 
system may be constructed by splitting a BECs, which 
is already coupled to a cavity mode, into three weakly 
linked condensates in many ways. In the large-detuning 
limit and in the rotating frame at the pump frequency, 
the Hamiltonian for the condensate system can be writ- 
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Figure 1: The sketch of cavity-mediated triple- well system 
in which (a) is the symmetry case (A = 0) and (b) is the 
asymmetry case (A < 0). A is the zero-point energy in each 
well and the BECs will be trapped in the wells. The cavity is 
driven by a pump field with an amplitude r\ and its loss rate 
is ft. 



ten as (h=l throughout), 

H = H a + Hf + Hi nt , 



(2.1) 



H a is the condensates Hamiltonian in the three-mode 
approximation, 

H a = \ib\h - b\b s ) - v(b\b 2 + b% + b\b 2 + b%) 

+ ic(6 t 1 6 t 1 6 1 6 1 + b\b\b 2 b 2 + btblbsbs), (2.2) 

Where b\(bi)(i = 1,2,3) creates (annihilates) an atom 
in its internal ground state in the ith traps, v is the tun- 
neling matrix element between the nearest modes, while c 
denotes the repulsive interaction strength between a pair 
of atoms in the same mode. A is the zero-energy of the 
wells. Hf is the external single-mode field Hamiltonian, 



Hf = ui c a*a + i](t)c 



rj*{t)e 11 



(2.3) 



where uj c and u> p are the cavity-mode frequency and the 
pump frequency, respectively. rj(t) is the amplitude of 
the pump field and varies slowly, i.e., \r]/r]\ <C uj p . In the 
limit of large detuning [U and the weak pump field, the 
atom-cavity field interaction is of dispersive nature, in 
this sense the upper level of the atoms can be adiabat- 
ically eliminated. Under the three-mode approximation 
for the BECs and dropping the coupling terms between 
the cavity mode and the atomic tunneling, the interac- 
tion Hamiltonian between BECs and the cavity mode can 
be obtained, 



Hi, 



UoctaiJtblb! + J 2 b\b 2 + J 3 blb 3 ). (2.4) 



where Uq = go 2 /{lo c — ui a ) characterizes the strength of 
the light shift per photon that an atom may experience, 
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with go being the atom-cavity mode coupling constant. 
The parameter J^(i = 1,2,3) denotes the overlap be- 
tween the atomic mode and the cavity mode. From Eq. 
(|2.4p , one may immediately find that the frequency cavity 
mode is renormalized due to the atom-field interaction, 
meanwhile the atom ensemble feels a induced potential. 
This atom-field nonlinearity will induce novel phenomena 
in the triple well system. According to the Heisenberg's 
equation, a dimensionlcss Schrodingcr equation of this 
cavity-mediated triple-well system can be written as 

161 = -vb 2 + A&i + cb\bibi + JitWa&i (2.5) 

162 = —vb\ + v b 3 + + J2Uoa [ ab2 
it>3 = — vbj, — A63 + 0636363 + J 3 XJ§o^ab 3 

ia = [uj c + UoiJiblb! + J 2 b\b 2 + J 3 b\b 3 )]a 

-iKa + r](t)e-' lu *\ (2.6) 

the cavity loss K is induced phenomcnologically, which 
is the dominant dissipation process since spontaneous 
emission is suppressed under the large detuning approx- 
imation. Eq.(2.5) shows that, the dispersive interaction 
between cavity photons and the condensate atoms intro- 
duces an effective light shift to the states 6,(i = 1,2,3). 
The characteristic values of k of the high-finesse opti- 
cal cavity, the order of 2ir x 10 6 Hz, is typically much 
larger than the spin oscillation frequency with the order 
of (2n x 10) — (27r x 10 2 ) Hz, This means that the cavity 
mode can follow the condensates wave functions adiabat- 
ically. Taking the standard mean-field approximation, we 
treat the operators bi(i = 1,2,3) and a as classical quan- 
tities, 61 = vOvTe 1 " 1 , 6 2 = VJhe i02 , 63 = ^/N 3 'e ie3 , and 
a ~ a. Here Ni(i = 1,2,3) are the numbers of atoms in 
the three wells, and 8i(i = 1,2,3) are their phases. The 
relative phases are 4>\ = 6\ — 62 and (f> 3 = 63 — 62- And 
the total number of atoms N = N t + N 2 + N 3 . Then the 
intracavity field amplitude can be derived from Eq.(2.6) 




Figure 2: The eigenenergy band structures of the condensates 
both for the uncoupled (a) and coupled (b-f ) cases. From (b) 
to (f), the reduced cavity detuning B is varied. The other 
parmeters are set as 5U A 2 /(2v) = 0.02, C = 0.07. 



reduced detuning, and the reduced loss rate, respectively. 
The photon number Eq. ()2.8[) can be written as 



HOI 5 



A{tf 



[m + n 3 + Bf + C 2 



(a) = a(t) = : 

-iK + [w P -cjc- Uo{JiN x + J2N2 + J 3 N 3 )] 

and the intracavity photon number is 



where rii=Ni/N(i = 1, 2, 3), and m 
the Eqs. (5) can be rewritten in tern 
the phase difference cf>i(i = 1,3) as 



(2.9) 

= 1. Then 
1, 3) and 



(aTa) = \a(t)\ 2 = 



W*)l s 



[A - J 2 NU + 5Uq(Ni + N 3 )] 



lil = 1 — n-1 — n 3 sin <f)\ , 



(2.10) 



where A = u p — uj c — J2NU0 and S = J2 — J\ (= J2 — J3) 
is the coupling difference between the two atomic modes 
to the same cavity mode. From Eq. (|2.8[) . one may eas- 
ily conclude that, the important feature of the cavity- 
induced effective light shift is that it is sensitive to the 
population distribution of the condensates in the triple 
well, which in turn will induce novel nonlinear phenom- 
ena in BECs. By introducing three dimensionless pa- 
rameters A(t) = rj(t)/[5NU ], B = A/[SNU ], and 
C = k/[SNUq] named the reduced pumping strength, the 



ih = V^Vl - ni - n 3 sin 03, 



(2.11) 



bi= -A - r[m - (1 - m - n 3 )} + I v - -j= cos0! 



1 "3 

/nT 



1 ( v /nTcos0i + v /nj"cos<fe) , Slip | a |2 (2 12) 

2 —713 2v < < > V ' 



03 = A - r[n 3 - (1 - m - n 3 )} + \ ^^f" 3 cos0 3 

+ 4^M 2 , (2.13) 



2 V"3 

1 {^/nTcostpi + ^/n^cos 4> 3 ) SUp | |2 
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r = 0, A = 2, the eigenenergy levels of the condensates (a) 
uncoupled and (b) coupled to a single-cavity mode given by 
Eq. (|2.14p and (c) the intra cavity photon number \a\ 2 versus 
the reduced cavity detuning B. 

Figure 3: For the parameters 5U A 2 /(2v) = 0.02, C = 0.07, 
r = 0, A = 0, the eigenenergy levels of the condensates (a) 
uncoupled and (b) coupled to a single-cavity mode given by 
Eq. (|2.14[) and (c) the intra cavity photon number \a\ 2 versus 
the reduced cavity detuning B. 



where the time is rescaled in units of the Rabi oscil- 
lation time l/(2i>), 2vt — > t. The dimensionless pa- 
rameter r = Ncj (2v) measures the interaction strength 
against the tunneling strength. From Eqs. ()2.12|) and 
(|2. 13|) . one can find that the cavity field influences the 
atomic dynamics in the same manner as the the zero- 
energy, i.e. A, of the triple wells system, which will lead 
to a redistribution of the atomic population among dif- 
ferent states hi based on the particle tunneling. Impor- 
tantly, this cavity-induced effective zero-energy is depen- 
dent on the atomic population distribution, as manifested 
by Eq. (|2.8[) . It is precisely this interdependence of 
the atomic dynamics and the intracavity photon num- 
ber that will result in some interesting nonlinear dy- 
namics of this cavity-mediated system, which is the fo- 
cus of this paper. Furthermore, a classical Hamiltonian 
H c (ni, n 3 , 0i, 03, t) can be constructed, in terms of two 
pairs of conjugate variables ri\ and </>i, 71,3 and (f> 3 , using 
rii = —jnf^, 4>i = TpS Such a counterpart mean- field 



Hamiltonian of the quantum Hamiltonian (|2.ip reads 

H c = \{n 3 - Tii) - \r[n\ + n§ + (1 - i%i - n 3 ) 2 } 
+V1 - n>i - n 3 (^/n^ cos 4>i + y^njj cosfo) 

-^F( ni ,n 3 ,t), (2.14) 

with 

A 2 (t) r ni + r2.3 + _B n 
F(m,n 3 , t) = — ¥- arctan[^ -± ]. (2.15) 

o o 

The first three terms in Eq. (|2.14p are the Hamilto- 
nian of a bare Josephson Hamiltonian as first derived 
in Ref. [3l|. They describe the energy cost due to the 
zero-energy, the phase twisting between the two conden- 
sates and the atom-atom repulsion, respectively The last 
term in Eq. (|2.14[) is the cavity-mediated interaction. In 
its nature, this term is similar to the potential felt as the 
atom passing trough the cavity field adiabatically j38j . 
which tilts the symmetric well. In this paper, we only 
concentrate on the case that the pump strength is a con- 
stant, so that the cavity- mediated triple well system is 
autonomous and the Hamiltonian (|2.14[) is conserved in 
time. 
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Figure 5: The eigenenergy levels as the parameter r is varied 
for 5U A 2 /(2v) = 0.02, B = -0.9, C = 0.07. 
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Figure 6: (Color online) The average of (m) (i = 1,2,3) 
of the BECs for the different interaction r, with initial value 
m(0) = l[(a) and (b)], ri2(0) = l[(c) and (d)] in the symmet- 
ric wells A = 0, respectively, (a) and (c) are for the uncoupled 
cases while (b) and (d) are for the coupled cases. The other 
parameters SU A 2 /(2v) = 0.02, B = -0.9 and C = 0.07 



III. THE ENERGY STRUCTURE 

In this section, we will study the eigenenergies of the 
condensates in the parameter regime that supports the 
existence and the conditions for loop structure emer- 
gence of the adiabatic lower and higher energy levels in 
the cavity-mediated triple well system and presents the 
unique features which make them different from those in 
the optical lattices and double well. In order to obtain 
the effects of the nonlinearity to be clearly identified, we 



chose c = firstly, which can be achieved via Feshbach 
resonances or is approximately valid in a sufficiently di- 
lute gas. The cases including the atomic collisions are 
also calculated in the last part of this section for com- 
parison with the former works [31| . The fixed point or 
minimum energy point of the classical Hamiltonian sys- 
tem (|2.14[) corresponds to the eigenstate of the quantum 
system. Solving the eigenvalue equation, one can derive 
the analytical expressions for these fixed points, how- 
ever, the analytical solutions are lengthy, it is difficult 
to interpret them, then we restore to the numerical and 
gra phic method to show the new physical effects in them 
J39|. The eigenenergies of the BECs uncoupled or cou- 
pled to the cavity field as a function of the zero-point 
energy bias A arc shown in Fig. [5J with Fig. 2(a) for 
the uncoupled cases, Figs.2(b)-2(f) for the coupled cases 
with the increasing dctunings. It can be found that the 
cavity-mediated BECs has a more complex and diverse 
energy structure. To be special, the coupled cases show 
the bistability behavior and the loop structure appears, 
which are not observed in the uncoupled counterpart, 
which implies that some new stationary points appear in 
the eigenenergies due to the cavity-mediated interaction. 
We can see that with the increase of the detunings, the 
lower energies E 3 and E 2 show bistability firstly, then the 
higher energy E\ does too. When the detunings are very 
large, the bistability behavior disappears for all the en- 
ergy levels. For the suitable detuings, the lower and the 
higher energies E% and E\ show loop structures, which 
shrinks in size and vanishes as the detuning is further 
increased. But the energy E 2 does not show loop struc- 
ture for the parameters considered in Fig. ([2]). Actually, 
numerical calculation for a wide parameters shows that 
the energy E 2 does not has a loop structure. Seen from 
Fig. 2(b)-2(f), it can be found that the emergence condi- 
tion for the bistability are different for different energies. 
i?3 only shows bistability phenomena at near A = as 
the detunings is increased, and loop structure disappears 
for large detunings. E 2 show bistability behavior at near 
A = when the detuning is relatively small, while for 
the large detunings, the bistability appears for large A. 
As for Ei, bistability begin to appear at the regions of 
larger A for the large detunings; With the increase of the 
large detuning, the bistability regimes move towards to 
the directions of the small zero-energy A. 

Looped structures have been discussed before for BECs 
in the ordinary (noncavity) f40M42| , the cavity- mediated 



optical lattices [21[ and the bare tripe well [31|, respec- 
tively, which is important due to the possible wide impli- 
cations for various systems in condensed matter physics 
and nuclear physics 



43] 



But the loop structure in Fig. 
([2]) is qualitatively different from the ones, resulting from 
the interatomic interactions nonlinearity obtained in Ref. 
j3l| . Here the bistability comes from the cavity- induced 
nonlinearity. It can be understood as follows: Because 
of the sensitive dependence of the BEC dynamics on the 
cavity induced potential. The cavity-induced phase shift 
can effectively result in population redistribution among 
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Figure 7: (Color online) The average population of (a), (c), (e) 
(uncoupled case) and (b), (d), (f) (coupled case) (rij) with the 
initial value m = 0.5, 712 = 0.5 and 113 — in the symmetric 
wells A = 0. The other parameters of the cavity are chosen 
as SUoA 2 /{2v) = 0.02, B = -0.9 and C = 0.07. 



Figure 8: (Color online) The average population of (a), (c), 
(e) (uncoupled case) and (b), (d), (f) (coupled case) (m) with 
the initial value ru = 0.5, 713 = 0.5 and 712 = in the sym- 
metric wells A = 0. The other parameters of the cavity are 
chosen as 6U A 2 /(2v) = 0.02, B = -0.9 and C = 0.07. 



the three wells even at the single-photon level. The BEC 
atoms in turn collectively act as a dispersive medium 
shifting the cavity resonance. Bistability results from this 
nonlinear feedback between the photons and the atoms 
in different wells. With the increase of the detunings, the 
cavity induced potential begin to enhance the triple-well 
asymmetry felt by the atoms. As a result, the atoms 
tend to stay in the deep well. This will give rise to a 
larger population imbalance and shift the cavity reso- 
nance, and will again lead to the increase of the probe 
mode intensity, as can be seen from Eq. (7) by setting 
detunings < —Uo(Ji,Ni). 

The trend of the energy structures and intra-cavity 
photons with respect to the reduced pump-cavity detun- 
ings are shown in Fig. [3] We focus on the condition 
A = and A = 2, respectively. Beginning from A = 
case, it can be found that, the level structure is simi- 
lar to the uncoupled counterpart for the small detuning 
(B < B\ = —0.994). While with an increasing detun- 
ings (£>i < B < B 2 — —0.849), the bistable curve ap- 
pears in the lower energy level, which corresponds to the 
loop structure of the lower level E3 in Fig. Htc). In the 
bistable regions, there are two minima and one maxi- 
mum that correspond to three different motional modes 



of the atoms in the third energy. As a consequence of the 
bistable curves or the loops structure, the existence of the 
multiple minima value at the fixed reduced detuning B is 
a omnipresent sign of hysteresis [H| . When the detuning 
is larger and exceed the critical value B3, the bistable 
curve emerges in the upper level, which corresponds to 
the loop structure of the upper level E\ in Fig. [He). The 
bistability can also be observed in the same ranges of de- 
tunings in the mean intracavity photon number |a| 2 , as 
shown in Fig. [3jc). The results for A = 2 cases arc shown 
in Fig. |H with the increase of the detunings, the energy 
E2 also shows bistability behavior between B\ = —0.96 
and B 3 = -0.3. 

Figures [U [21 and 2] focus on the non-interaction cases. 
Now we turn to the interacting BECs case. Figure [S] de- 
picts the eigenenergy versus the zero-point energy bias A 
but covering different interaction strengths r. As we con- 
sider two types of the nonlinearity simultaneously, i.e., 
the atomic interaction and cavity-mediated nonlinearity. 
the energy structure is more complex than the one ob- 
tained for only one type of nonlinearity mentioned. It is 
shown that the structure is similar to the noninteraction 
case (r = 0) when the parameter r is weak (r = 1). A 
loop structure appears at the middle level E2 shown in 
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Figure 9: (Color online) The average population of (a), (c), 
(e) (uncoupled case) and (b), (d), (f) (coupled case) (n») with 
the initial value ni = 1/3, ni = 1/3, 713 = 1/3 for different 
relative phases in the symmetric wells A = 0. The other 
parameters of the cavity are chosen as SUoA 2 / (2v) = 0.02, 
B = -0.9 and C = 0.07. 



Fig. [5fb) as r is increased (r~l). As r increases to 1.5 
and 2.5, respectively, the loop structures have more com- 
plex shapes in the upper level E\ and the middle level E2 
as shown in Figs. [5jc) and[5jd). Compared Fig. [5] with 
Fig. [3J one can find that, this kind of complicated loop 
structure occurred in the upper and middle levels comes 
from the nonlinearity caused by the atom-atom interac- 
Furthermore, the loop structure of the lower 



tion 31 



level E% does not change with the increase of the inter- 
atomic interaction as shown in Fig(5] Thus one can see 
that this loop structure in energy E3 is mainly caused by 
the nonlinearity generated from the cavity back-action 
rather than from the direct interatomic interactions. 



IV. TUNNELING TO LOCALIZATION, A = 

In this section, we will perform numerical simulation 
of localization or self trapping of a BEC in this cavity- 
mediated triple well system using the solution of the Eqs. 
()2.10j) - ()2.13[) . In general, the localization or self trapping 
is caused by the nonlinear interaction between atoms, 
then the self trapping of a repulsive BECs system in a 



Figure 10: (Color online) The average of (m) (i — 1,2,3) 
of the BECs for the different interaction r in the tilted wells 
A = — 1, with initial value ni(0) = l[(a) and (b)], 712(0) = 
l[(c) and (d)], and 713(0) = l[(e) and (f)] respectively, (a), 
(c) and (e) are for the uncoupled cases while (b), (d) and (f) 
are for the coupled cases. The parameters of the cavity are 
SU A 2 /(2v) = 0.02, B = -0.9 and C = 0.07. 



bare optical lattice, double well and triple well potential 
has been investigated firstly both for the dynamical and 
stationary cases [Hj]. When considering the coupling to 
a cavity, the cavity-mediated atoms-field interaction will 
provide a novel nonlinearity different from the atoms- 
atoms nonlinear interaction, then the research shows the 
self-trapping behavior can be modified drastically by the 
strong condensate-field interplay @, [H, [I?} in a double 
well potential. As mentioned before, the cavity-mediated 
triple- well system induce a effective potential and tilt the 
symmetric well as a result. The more asymmetric the 
well becomes, the more complex the dynamical behavior 
appears. Then it is interesting to investigate the local- 
ization and the self-trapping dynamics in this system. 
We will calculate the average population of each well for 
the different initial values with a different parameter r 
in the following. As a benchmark, we firstly discuss the 
symmetry cases 7 = 0. 
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Figure 11: (Color online) The average population of (a), (c), 
(e) (uncoupled case) and (b), (d), (f) (coupled case) (nt) with 
the initial value ni = 0.5, ri2 = 0.5 and 713 = in the tilted 
wells A = — 1. The other parameters of the cavity are chosen 
as SU A 2 /{2v) = 0.02, B = -0.9 and C = 0.07. 



Figure 12: (Color online) The average population of (a), (c), 
(e) (uncoupled case) and (b), (d), (f)(coupled case) (n*) with 
the initial value n\ = 0.5, nz = 0.5 and n% = in the tilted 
wells A = — 1. The other parameters of the cavity are chosen 
as 5U A 2 /{2v) = 0.02, B = -0.9 and C = 0.07. 



A. localization in one well 

Comparing to a bare BECs in the triple-well, here the 
BECs in the cavity-mediated triple well systems couples 
to the cavity mode, the atom-cavity interplay terms, the 
last term of the Hamiltonian (|2.14[) , will not only change 
the critical value r c , but also alter the process of transi- 
tion to localization. To illustrate this point, we calculate 
the mean value (m) (i — 1,2,3) of the BECs uncoupled 
or coupled to the cavity mode for different interaction 
r with the initial conditions rii(0) = 1 and 712(0) = 1, 
respectively. The results are presented in Fig. [51 When 
the BECs are uploaded initially in the left well, shown in 
FiglUb), it can be found that the average populations of 
each well are more chaotic than the ones of the uncoupled 
case (Fig. [H^a)), which is due to the reasons that the non- 
linearity of the effective atom-cavity induced potential 
in Eq. (|2.14[) tilts the wells and changes the stationary 
points of the condensates and makes some new motional 
modes of the atomic oscillation appear, correspondingly 
the critical value of the transition in the coupled case is 
larger than the uncoupled one. However, when the BECs 
are uploaded in the middle well initially, shown in Fig. 
in^c) and^d), the cavity-mediated nonlinearity does not 



have a pronounced effects on the transitions. It can be ex- 
plained as follows: as mentioned before, the atom-cavity 
induced interaction is equivalent to the zero-energy A, 
which affect the dynamics by tilting the symmetric well. 
When one put atoms in the middle well, the asymme- 
try of the wells 1 and 3 does not affect the dynamics 
of the BECs in well 2 apparently With the increase of 
the interaction, the oscillation will be compressed and 
the BECS atoms will be trapped in the middle well com- 
pletely. Comparing the coupled cases with uncoupled 
cases, one can also find that the values of {n\) and (712) 
are different for r = 0, which is implied that the BECs 
atoms in the wells are rearranged when the condensates 
couple to the cavity. 



B. localization in two- and three- wells 

This part will focus on the dynamics of the BECs 
atoms trapped initially in two- and three-wells with dif- 
ferent initial relative phase. We calculate the average 
value (m) (i = 1, 2, 3) as a function of r for different cases. 
(rii) versus r with the initial value n\ = 0.5, n-i = 0.5, 
and 713 = (the BECs atoms are initially uploaded in 
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Figure 13: (Color online) The average population of (a), (c), 
(e) (uncoupled case) and (b), (d), (f) (coupled case) (n») with 
the initial value ni = 1/3, ni = 1/3, 713 = 1/3 for different 
relative phases in the tilted wells A = — 1. The other parame- 
ters of the cavity are chosen as 5UqA 2 / (2v) = 0.02, B = —0.9 
and C = 0.07. 



wells 1 and 2) is shown in Figs. [7] For comparison, the 
uncoupled cases are also presented. As BECs are up- 
loaded in two wells initially, the interference should be 
considered, and their initial phase will be important. In 
general, different initial phase will lead to different behav- 
iors. From Fig. one can find that the coupling to the 
cavity will generate distinct effects for different relative 
initial phase, when the relative phase is zero, the atom- 
cavity coupling will make the smooth oscillation process 
become chaotic process for the small interaction r, as a 
result the coupling to the cavity will enlarge the chaotic 
region. While for the relative phase 7r/2 cases, this atom- 
cavity nonlinearity will decrease the chaotic region. With 
the increase of the interaction, the atoms will be trapped 
in wells 1 and 2 for the relative phases 7r/2 and zero, re- 
spectively. When the relative phase is set to 7r, compared 
Fig- with Fig. UJe), one can find that the atom will 
be trapped in wells 1 and 2 from small interaction r, but 
the atoms in each well oscillate quickly and the atoms 
can not be trapped easily in each well. Furthermore, the 
cavity interaction will change the population in wells 1 
and 2. In the uncoupled cases, the population (ni) is 
higher than 0.5 for a large interaction region, the cou- 



pling to the cavity will lead the population (rii) decrease 
to a lower value than 0.5 for a large interaction region as 
shown in FigJTJb). While for (n- 2 ), the coupling to the 
cavity will lead the population (712) to be higher than 0.5 
for a large interaction region on the other contrary. Com- 
paring FigJTJd) with Fig. [5] (d), one can also find that 
the BECs in the middle well are more chaotic as atoms 
trapped in two wells initially due to the tunneling and 
interference coming from the BECs in the well 1. The 
mean value (rii) versus r with the initial value n\ = 0.5, 
?i2 = 0, and n 3 = 0.5 (the BECs atoms are initially up- 
loaded in wells 1 and 2) is shown in Fig. [5J It can be 
found that the atom-cavity coupling results into differ- 
ent effects for the different relative phase as mentioned 
in [7J After passing a wide oscillation region, the oscilla- 
tion and the chaos behavior are compressed into a narrow 
region and the atoms will be trapped equally for a large 
interaction r in the two traps eventually where they are 
initially uploaded for three cases. From Fig(5Jb), we can 
find that the noninteracting atoms are mostly trapped in 
well 1 due to the coupling to the cavity for the relative 
phase 7r/2, but with the increase of the interaction, the 
population decrease quickly and enter a chaotic process. 
When the relative phase is chosen as zero, the cavity- 
mcdiatcd effects only enlarge the oscillation region, but 
for the other two relative phases cases, this nonlinear- 
ity not only enlarge the oscillation region but also make 
more chaotic behavior. Then we can change the relative 
phases to control the dynamics of the BECs in this cavity- 
mediated triple well. The above are scheduled to study 
the localization of BECs in two wells. As the triple well 
system is a protype of optical lattices helping us provide 
a bottom-up understanding of mechanisms operating in 
the infinite optical lattice. Here, we will turn to the case 
that the BECs atoms are trapped in three wells simulta- 
neously with initial n\ = 1/3, 712 = 1/3, and 713 = 1/3. 
We calculate the average value (ni) (i = 1,2,3) as a 
function of r for different relative phases in Fig. [S] The 
uncoupled cases are also shown for comparison. It can 
be easily found that the dynamics is also sensitively re- 
lated to the initial phase but the phase does not show a 
pronounced effect as it does in the two wells cases. The 
coupling to the cavity also enlarges the oscillation re- 
gion and make more chaotic for the cases of the relative 
phase 7r/2 and tt. While the relative phase is zero, the 
atom-cavity interplay may be neglected for its very weak 
influence on the dynamics, the BECs atoms show same 
process moving smoothly to a robust trapped behavior 
from the a small interaction r both for the coupled and 
uncoupled cases. 



V. TUNNELING TO LOCALIZATION, A / 

In the above, we have studied the dynamics of the 
cavity-mediated BECs in the symmetric triple- well cases 
(A = 0). Though the cavity- mediated potential makes 
the triple-well tilted as the zero-energy A does, the ac- 
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tual experiments values, SUqA 2 / (2v) — 0.02 << 1 only 
leads the triple well to tilt slightly. How does the cavity- 
mediated potential affect the dynamics of the BECs for 
the asymmetric triple- well cases, and what will happen 
for a cavity-mediated tilted triple-well system? In order 
to answer these interesting questions, this section will 
extended to the general case with A = — 1. In the follow- 
ing, the self trapping or localization in one-, two-, and 
three-wells are investigated, respectively. 



A. localization in one well 

We explore the dynamics of of the BECs by depict- 
ing the mean populations for the uncoupled and cou- 
pled cases versus the parameter r with A = — 1 and 
the initial conditions are rii(0) = 1, 712(0) = 1, and 
n 3 (0) = 1 , respectively. The numerical results are shown 
in Figs. fT0la)-fT0Ff). In Fig fTUT b). with an increase in 
the interaction r, the averaged population (ni) decreases 
smoothly at the beginning, passes a large oscillating and 
chaotic areas till r = 7.61, and then tends to unite. 
Compared to the uncoupled case (Fig. fTCT a)). it can 
be found that the transition area and the critical value 
r are larger for the coupled case due to the reason that 
the atom-cavity interplay enhances the oscillation of the 
BECs atoms of each wells. Thus a larger interaction r is 
used to trap the condensates. Compared FiglTOTb) with 
Figin^b), one can also find that, when the interaction 
is small, the coupling to the cavity will help the BECs 
atoms nearly be trapped in well 1 for the tilted case. 
In Fig. 110( d). the cavity-mediated interaction will make 
the averaged population (ni) show a turbulent behavior 
from a small interaction r, passes a wide chaos region, 
and then jumps up and tends to unite. In Fig. 110( f). ini- 
tially, upload all the atoms in the right-hand well, with 
an increase in r, the average population (713) undergo a 
smooth process, then tend to unite. In comparison to 
Fig. not e), it is shown that the coupling to the cavity 
will lead to a larger critical value of the interaction r. 



B. localization in two- and three- wells 

Except the wells are tilted, the parameters and the 
initial conditions in FigJTT] are as same as the ones in 
FigEl Compared with these two figures, we can find that 
besides the phase sensitivity and the enlarging chaotic 
region, the difference is, the condensates atoms for the 
relative phase it will be trapped too for the tilted wells 
cases. Furthermore, the mean value of (ni) is always 
small than (712) the due to the tilted wells, i.e., the atom- 
cavity induced tilt potential is weak than the tilt of the 
well, and then it can not make a inversion as it does in the 
symmetric well cases. The mean value (rii) for different 
initial phases with the initial value n\ = 0.5, n 2 = 0, 
n 3 = 0.5 for the tilted well cases is shown in FigTJ"2"1 
Compared to Fig. [51 we may find two characteristics 



due to the tilted well, one is the the smooth motion of 
atoms for the relative phase zero become more chaotic 
due to the common effects of the tilted well and the atom- 
cavity induced potential. The other is the trapped mean 
population of (n\) is large than (rii) due to the tilt. For 
comparison with Fig. [SJ the BECs atoms trapped in the 
tilted three wells simultaneously with initial n\ = 1/3, 
7i-2 = 1/3, and 773 = 1/3 are also investigated shown in 
FigHSl One obvious phenomena is for the relative phase 
zero case, when the well is tilted only, the mean values 
of rii is still smooth. When the tilted tripe well couples 
with the cavity, the atom-cavity interaction will destroy 
this Smooth behavior and make the dynamics become 
chaotic for small interaction r. 



VI. SUMMARY 

In this paper, we have studied the bistable curves and 
the tunneling dynamics of BECs in the cavity-mediated 
triple-well system within the mean field approach. We 
find that the atom-cavity nonlinearity will make the 
bistability occur both in the number of the photons of 
the cavity and in the energy structure of the BECs. Fur- 
thermore, the loops appear for the upper and lower en- 
ergy levels with the pump-cavity detuning increasing. We 
also demonstrated that the transition to localization of a 
BECs can be efficiently controlled by the optical cavity. 
To understand the effects resulting form the cavity cou- 
pling, we have performed a detailed numerical simulation 
both for the symmetric and asymmetric cases. 

For the cases of BECs atoms uploaded in well 1 ini- 
tially, the atom-cavity field nonlinearity enlarges the os- 
cillation and chaotic regions both for the symmetric and 
asymmetric cases, then a large interaction is needed to 
trap the atoms. When the BECs atoms are uploaded in 
well 2 initially, this nonlinearity only have a pronounced 
effects for the asymmetric case, which makes more chaos. 

For the cases of BECs atoms uploaded in two wells 
initially, the cavity have different effects for different rel- 
ative phases. With the increase of the interaction, the 
atoms will be trapped in one or two wells where they 
initially uploaded. Besides this, the coupling to cavity 
can lead a mean population inversion for certain relative 
phase if the BECs atoms are initially uploaded in wells 1 
and 2. 

For the cases of the BECs atoms uploaded in three 
wells initially, there is a robust trapped phenomena for 
the relative phase zero in the symmetric wells even for 
the coupled cases. The theoretical discussion here will 
be helpful in controlling of BECs in real experiments. 
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